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ABSTRACT 

How do magnetohydrodynamic waves travel from the fully ionized corona, 
into and through the underlying partially ionized chromosphere, and what are 
the consequences for solar flares? To address these questions, we have developed 
a 2-fluid model (of plasma and neutrals) and used it to perform ID simulations 
of Alfven waves in a solar atmosphere with realistic density and temperature 
structure. Studies of a range of solar features (faculae, plage, penumbra and 
umbra) show that energy transmission from corona to chromosphere can exceed 
20% of incident energy for wave periods of one second or less. Damping of waves 
in the chromosphere depends strongly on wave frequency: waves with periods 10 
seconds or longer pass through the chromosphere with relatively little damping, 
however, for periods of 1 second or less, a substantial fraction (37%-100%) of 
wave energy entering the chromosphere is damped by ion-neutral friction in the 
mid and upper chromosphere, with electron resistivity playing some role in the 
lower chromosphere and in umbras. We therefore conclude that Alfvenic waves 
with periods of a few seconds or less are capable of heating the chromosphere 
during solar flares, and speculate that they could also contribute to electron 
acceleration or exciting sunquakes. 

Subject headings: magnetohydrodynamics (MHD) — plasmas — Sun: chromosphere 
— Sun: corona — Sun: flares — waves 



- 3 - 



Introduction 



During solar flares, up to 10^^ J of energy is released in the corona by reconflguration 
of magnetic flelds and a signiflcant fraction of this energy is transported ~ 10^ km to the 
chromosphere to power the UV and white hght emissions that make up the majority of flare 
electromagnetic radiation. This energy transport is widely ascribed to beams of electrons 
that are assumed to be accelerated in the corona in the vicinity of the magnetic energy 
release. Such models are appealing because electron beams impacting on the chromosphere 
can explain t he hard x-ray (HXR) source s associated with flares, using the coUisional thick 



target model (IBrown 



1971 



Hudson 



1972|). 



Electron-beam-only models face a numb e r of theoretical difficu lt ies w h ich have been 
descri bed elsewhere, e.g. 



Brown et al. 



fll990h : Fletcher and Hudson! (120081 ): 



Brown et al. 



( 120091 ). For example, coronal acceleration models require large numbers of particles when 
compared to typical coronal densities, or else a powerful resupply mechanism which is 
opposed by microinstabilities. Further difficulties include the small area of HXR footpoint 
sources determined by RHESSI observations, whic h imply electron beam densities in the 



corona high enough for the beams to be unstable (iBrown and Melrose 



1977 



Krucker et al 



20111 ). It is therefore worthwhile to consider other mechanisms by which at least some 
portion of flare energy could be transported to the chromosphere and deposited there. 

An alternative transport mechanism that has lately received renewed attention is 
magnetohydrodynamic (MHD) waves. Recent advances in observations have led to a 



consensus that the so 



ar atmosphere is 



pervaded by MHD wayes, which are now observed 



almost routinely (e.g. 


Ofman et al. 


1997 




Deforest and Gurman 


1998; 


Thompson et al. 


1999; 


Rob ; 


recht et al. 


1999 


. Aschwanden et al. 


1999 




Nakariakov et al. 


1999 




De Pontieu et al. 




2007 


Okamoto et al. 


2007; 


Tomczvk et al 


. 2007 


Jess et al. 


2009 


Mcintosh et al. 


2011 


, and 



references therein). Against this background, it has been proposed that a substantial part of 
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the energy released by a flare is in the for m of MHD waves, providing an alternative mode 



20081 ). possibly operating 



of energy transport to the chromosphere (IFletcher and Hudson 
alongside electron beams. 

There are a number of theoretical and observational arguments for considering 
MHD waves in this context, besides their ubiquitous nature. Theoretical studies of 
time-dependent magnetic reconnection under coronal conditions show that Poynting flux 
i s the dominant downwards energy flux beyond the diffusion reg i on and reconnectio n jets 



(ILongcope and Priest 



2007 



Birn et al. 



2009 



Kigure et al. 



20101). 



Birn et al. 



(120091) also 



demonstrated that Poynting flux can be particularly intense in tightly focused spots that 
resemble the brightest chromospheric emission sources, the spots being associated with 
fragmentation of a 3D coronal current sheet. Furthermore, the energy flux carried by 
MHD waves in active regions can be comparable to the flare energy flux for reasonable 
perturbations because Alfven speeds can be large (up to a few times 0.1c) according to 



coronal field strengths deduced from gyros ynchrotron emission (IBrosius and White 



obtained by magnetic field extrapolation fiRegnier et al. 



from observations were discussed in detail by 



20061) or 



20081) (estimates of Alfven speeds 



Fletcher and HudsonI ( 120081 )). Magnetic field 



strengths and densities vary throughout active regions, and for most active regions the figure 
of a few times 0.1c refers to places where the Alfven speed is particularly high, however, 
these are also the parts of an active region where flare energy transport by Poynting flux is 
likely to be of most interest. 

That the magnetic field can be strongly perturbed during a solar flare is not in question: 
there are many observations of changes to both the line-of-sight and horizontal magnetic 



field at the photosphere a t the time of t 



2005 



Wang and Liu 



2010 



Su et al. 



l e flare impu 



2011 



Wang et al. 



sive p hase (e.g. 



Sudol and Harvey 



20121). A broader question is what 



effect such perturbations, which are presumed to be driven by magnetic rearrangement in 
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the corona, could have on the flare atmosphere. 



MHD waves have l ong been contemp 
the chromosphere, with 



ated to play a role in the energy balance of 
PiddingtonI (jl956l ) realizing that ion-neutral collisions in the 
chromosphere produce a wave damping that is many orders of magnitude greater than that 
produced by either viscosity or classical resistivity. This was first used to explain heating of 
the chromosphere by upgoing waves produced in the up per convection zone a n d pho tosphere 



(lOsterbrock 



196ll ). but later applied to solar flares by 



Emslie and Sturrockl (jl982[ ). who 



studied how Alfven waves might contribute to the heating of the tem perature minirnum 



region (TMR) that is inferred from optical and UV flare obs ervations (IMachado et al. 



Cook 



1979 



Cook and Bruecknei 



MHD, WKB treatment. 



1979 



Metcal 



Emslie and Sturrockl ( 



• et al. 



1978 



1990al ). Based on their single-fluid 



19821 ) concluded that Alfven wave damping 



is an energetically viable mechanism for TMR heating. 



The o rigin of tempera t ure ra inim um heating an d 



explained. 



Machado et al. 



(119781) and 



of the candidate mechanisms, and 



Metcalf et al. 



optica l flare emission is still not fully 



Metcalf et al.l (|l990bt) provide excellent summaries 
f l990b[ ) then rules most of them out 



based on ob servations of 5 flares. Dire c t exci tatio n of the TMR by electron beams , as 



proposed by Aboudarham and Henouxl (119861 ) and lAboudarham and Henouxl (Il987l ). or 



by proton beams, is ruled out on the basis of inconsistency with the observed hard X-ray 
and 7-ray signatures. TMR heating by photoelectrons generated by EUV or soft X-ray 
photoionization of the TMR by an overlying coro nal or transition-regio n source is also ruled 



Metcalf et al. 



(Il990b( ) as chromospheric 



out. The most likely mechanism is identified by 
backwarming, in which enhanced chromospheric continuum radiation produced by the flare 
is directly absorbed (by the H~ ion) in the TMR, heating it and producing enhanced optical 
emission from the TMR and from the photosphere below (itself heated by the increased 
TMR emission). This mechanism is supported also by numerical simulations based on 



- 6 - 



energy deposition and excitation by an electron beam from the corona (e.g. ICheng et al. 



20101 ) and by observational compa risons with calcu lations involving electrons accelerated 



more locally to the radiation site ( iDing et al 



19991 ). However the role played by Alfven 



waves in chromospheric heating has not been explored to a comparable degree, so as part 
of this study we revisit the question of TMR heating by this mechanism. 



by 



It should b e note d that the wave dissipation model for TMR he ating was dismissed 



Cheng et al 



Emslie and Sturrock 



( 119881 ) on the basis that the time delay predicted by 
(Il982l ) between the appearance of flare hard X-rays and the TMR heating, due to the wave 
travel time from corona or upper chromosphere to the TMR, was too long compared to 



observations. On the other hand, 



Metcalf et al. 



(Il990bl ) ruled out the model on the basis 



that the predicted delay was a factor 2 too short to explain observations. The time delays 
depend on the local Alfven speed in the chromosphere, which may be expected to vary 
strongly from place to place in an active region, and is dependent on the vertical variation 
of the chromospheric density and field strength, neither particularly well constrained by 
observations or models. 

Typically, TMR temperature incr eases of between 100 and 3 K compared to typ ical 



quiet Sun values have been reported (IMachado and Linsky^ il975i : iMetcalf et al. 



mm) 



located at a few hundred kilometers above the photosphere. The total energy flux required 

:he order of 10^ J m~^ s~^ (deduced by 



Watanabe et al. 



(|2010h bv an e stima te 



Fletcher et al. 



fl2007h 



to power the optical emission from a large flare is of ' 
dividing the radiated power in the flare reported by 
of the flare's optical footpoint area from their Figure 2, while from 
we can set more modest (lower) limits of around 10^ J m~^ s~^, based on the measured 
UV footpoint area (which is, however, typically larger than the optical area). These are 
therefore the typical values that must be provided by the wave Poynting flux entering and 
dissipated in the chromosphere. 
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With solar flares as motivation, this paper examines Alfvenic wave propagation from 
the fully ionized corona into and through the underlying partially ionized chromosphere. 
Particular attention is given to the transmission of energy from the corona to the underlying 
atmosphere and to damping of waves in the chromosphere through ion-neutral friction and 
neutral-enhanced electron resistivity. These are crucial steps towards detailed understanding 
of the role MHD waves play in flares. 

Recently, the propagation, damping and energy flux of Al fvenic waves in par- 



tially i onized plasma has b e en examine d in a solar c o ntext 



(2003 


); 


Khodachen 


<.o et al. 


Pandev and Wardle 


( 


2008); 



(2004, 



20061 



Leake et al 



Vranjes et al. 



mm; 



(20051 



b y iKumar and Roberts 



Pandev et al 



Zaqarashvili et al. 



(120081); 



(120111 ). confirming 



earlier results that ion-neutral collisions are responsible for significant dissipation in the 



chromosphere. Effects of partial ionizatio r i have also 



waves m prommences 
form a tion of spicules 



2002 



James et al. 



Forteza et al 



Haerendel 



20031). 



2007 



1992 



Soler et al. 



jeen considered for damping of 



2009 



2010 



de Pontieu and Haerendel 



Carbonell et al 



1998 



2010h. 



James and Erdelyi 



(IKhomenko and CoUados 



chromospheric heating by dissipation of magnetic flux tubes 



between ions and neutrals ( Krasnoselskikh et al. 



20121). and generation of ele ctric currents by relative flows 



2010f ). 



neutr als may be important for ma gnetic flux emergence (ILeake and Arbei 



20071) and rnagne tic reconnection (jZweibel 



Leake et al 



20121 ) 



1989 



Additionally, the presence of 



Sakai and Smith 



20091 : 



2006 



Arber 



Zweibel et al 



et al. 



2011 



Transmission of Alfvenic waves from the corona to the chromosp 
subseq uent damping have previously been investigated analytically by 



lere and their 



De Pontieu et al. 



( I2OOII ). who evaluated damping times at different heights in model chromospheres and 
computed energy transmission coefficients for a continuous Alfven speed profile that is 
uniform in the photosphere and corona, and exponentially growing with height in the 
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chromosphere, using a constant scale height. However, they did not include a step in Alfven 
speed at the transition region, and therefore overestimated transmission from the corona to 
the chromosphere. 

We revisit the subject using a two-fluid model (plasma and neutrals) that avoids 
adoption of the WKB approximation or assumptions of Cowling conductance. We also go 
beyond previous investigations by solving the complete time evolution of Alfvenic wave 
pulses launched from the corona, using realistic de nsity and temp e ratur e structure from a 



Fontenla et al. 



( I2OO9I ). This approach 



selection of semi-empirical atmospheric models by 
allows the most detailed and realistic investigation of the problem to date. Although this 
paper focuses on downgoing waves, the models and methods described can equally be used 
in understanding upgoing waves and may therefore be valuable to study of the solar wind, 
coronal and chromospheric heating, and spicules. 



2. Energy Transmission at the Transition Region 

Before beginning detailed investigation, it is useful to estimate energy transmission 
across the solar transition region for a simplified model where an Alfven wave encounters a 
step in Alfven speed, corresponding to the transition between a uniform corona and uniform 
chromosphere. Although highly simplified, this model produces some useful insights. 

In the corona, the wave speed is 

Bn 



cor 

where Bo is the background magnetic field strength and pcor is the mass density of the 
coronal plasma. For the chromosphere, we assume that neutral particles and ions are 
perfectly coupled by collisions, making the wave speed 

^y^J'OPchr 
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where pchr is the total mass density (ions and neutrals). 

Solving the linear momentum and ideal induction equations, subject to requirements 
that an infinitesimal region around the transition region has finite acceleration and is 
neither a source nor sink of wave energy, the energy transmission coefficient (fraction of 
incident wave energy that is transmitted) is obtained as 



Emslie and Sturrocki (jl982l ) in introducing a 



Here, we have followed the example of 
dimensionless parameter, 0, defined as 

I / C-chr \ Pcor ^cor / < \ 

= ( = = — - — , (4) 

V Ccor / Pchr ~r 

where the last equality holds for hydrogen-dominated plasma with ricor the coronal plasma 
number density, rig the chromospheric plasma number density, and Ua the chromospheric 
neutral number density. These results hold for waves incident from either the corona or 
chromosphere. 

It is convenient to express energy transmission in terms of the chromospheric and 
coronal temperatures. Assuming pressure balance at the transition region, 

'^^corTcor ('^a ~l~ '^^e)'-^chr- (5) 

Hence, 

/T , \ 

f{njna) (6) 



chr 

f 

cor 

where f{ne/na) = (0.5 + ne/na)/{l + ne/ua) is a monotonically increasing function, 
bounded between 0.5 (neutral chromosphere; rie/ria — )■ 0) and 1 (fully ionized chromosphere; 
^e/'^a oo). The parameter is typically of order 0.01 and can therefore be neglected on 
the denominator of (|3]). Thus, the energy transmission coefficient, T, lies in the range: 



2V2J^< T <4J^, (7) 



T V T 

-'- cor V -'■ cor 
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with T approaching the upper hmit if the upper chromosphere is fully ionized, and the 
lower limit if the upper chromosphere is dominated by neutrals. 

The limits on T are easily evaluated for typical solar conditions. For example, 
substituting T^hr = 10000 K and Tear = 2 MK into ^ gives the range 0.20 <T < 0.28. We 
therefore expect significant transmission of wave energy across the transition region. This is 
especially true if coronal waves are trapped (on closed field lines, within density structures 
or inside the chromosphere) and therefore incident on the transition region multiple times. 

The exact result is dependent on the temperatures chosen, however the dependence 
is lessened by the square root: even if the temperature ratio were increased fourfold, the 
energy transmission coefficient would still be between 0.10 and 0.14. If collisional coupling 
between chromospheric ions and neutrals were imperfect (which may be the case for high 
frequency waves) the primary effect would be to increase Cchr through reduced mass loading 
and hence increase transmission. We also note that, although the present analysis is founded 
on Alfven waves, it seems likely that similar conclusions will apply to other MHD waves 
(for example, fast or kink waves), something that should be confirmed by future work. 

We conclude from this short analysis that waves with an Alfvenic nature can carry 
significant energy across the transition region. This motivates further study, so we proceed 
to detailed modeling of Alfven waves in a solar atmosphere with realistic height structure. 

3. Modeling Alfven Waves in the Chromosphere 
3.1. Chromospheric Properties 



Semiempirical models for ch romospheric composit ion, densities and te mperature 
include the well-known models of IVernazza et al.l (119811 ) and, more recently, iFontenla et al, 
(120091 ). The latter have state-of-the-art physics and distinguish between different solar 
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features such as umbra, penumbra, plage and faculae. We have chosen non-flaring models 
because we are interested in how the chromosphere might get into the flaring state from a 
preflare condition. 



Figure [T] (top) and (middle) show particle numb er densities and 



Fontenla et al. 



temperature as 



(120091). Neutral 



functions of height for the plage model described by 
hydrogen is the most common particle and it is essentially in thermal hydrostatic 
equilibrium, although the model does consider parametrically-described non-gravitational 
forces. The neutral scale height is typically between 120 km (lower chromosphere) and 
230 km (upper chromosphere). Electrons and ions are also present, with the fraction of 
hydrogen that is ionized varying from 5 x 10^^ at the TMR to 1 at the transition region. 
At the TMR, Figure [1] (top) shows that electron number density is an order of magnitude 
greater than proton number density. The difference is made up by 'heavy' ions that 
substantially outnumber protons at these heights because hydrogen ionization is so low. 
Thus, magnesium, silicon and iron ions - the most abundant elements in the Sun with first 
ionization potentials lower than hydrogen's - should be included in chromospheric models 
that include the TMR. 

Magnetic field strength as a function of height in the chromosphere is not well 
established, due to the difficulty of making direct chromospheric measurements, but it is 
widely agreed that magnetic field strength in flux tubes should decrease with height. In the 
absence of observationally-derived models, we follow an established precedent and model 
the decay in magnetic field strength as 



Biz) = 5(0) 



P{0)J 



where P{z) = kTijia + 2ne) is the total gas pressure (iZweibel and Haberl Il983l ). For 



a plage model, we take B = 0.05 T at the photosphere, in keeping with magnetogram 
measurements, and B = 0.01 T at the top of the chromosphere. Combined with gas 



pressures from iFontenla et al.l (120091 ) . this gives the magnetic field profile shown in Figured] 
(bottom), for which a = 0.139. 



3.2. Coronal Properties 

Coronal properties partially determine wave transmission across the transition region, 
and including a coronal domain in our model also allows propagation of an incident wave 
to and from the transition region. Neither of these roles require realistic height dependence 
away from the transition region, hence we will model the corona as a uniform fully ionized 
plasma with properties appropriate to the low corona. Throughout this work, a coronal 
temperature of 1 MK is assumed, and we impose a condition of gas pressure balance at 
the transition region, from which coronal plasma density is deduced. The resulting coronal 
density varies between solar features from a few times 10^^ m~'^ to 10^^ m^^ (over sunspots), 
in reasonable agreement with accepted values for the low corona. 

We have also run additional simulations with lower coronal densities (assuming coronal 
temperatures of 2 MK and 4 MK). Halving coronal density reduces transmission to the 
chromosphere by a factor of \/2, in line with equation (^^, and does not alter our damping 
results. 

3.3. Governing Equations 

Section 13.11 shows that the chromosphere must be described as a partially ionized 
plasma with neutral and plasma components. One way to do this is using a multi-fiuid 
approach, which we outline here. 

The neutral component of the plasma may be modeled as a fiuid with the following 
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Neutral Hydrogen 

Electrons 

Protons 




0.6 0.8 1 
Height (Mm) 



Fig. 1. — Plage chrom ospheric model, (top) Number densities of neutral hydrogen, elec- 



trons and protons, from 



Fontenla et al. 



(120091 ). Where proton number densities depart from 



elect ron number densities, this is due to the presence of heavy ions, (middle) Temperature, 



from 



Fontenla et al. 



( 120091 ). (bottom) Magnetic field strength in a representative flux tube. 



assuming power law scaling with gas pressure. 
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momentum equation: 

Pn—^ = -^Pn + PnG + pnT^n.i {U-l - Un) + Pnl^n,e (^e " , (9) 

where Dg/Dt = d/dt + Ug ■ V is the convective derivative. Terms on the right-hand side 
represent, from left to right, pressure gradients, gravity and momentum exchange with 
charged particles via collisions. Here Ug^t are momentum transfer collision frequencies, 
generally a function of temperature and the number density of target particles, rit. These 
have the property 

PsJ^s,t = PtVt,s- (10) 

Momentum transfer takes the form of a friction, as shown in equation ([9]), provided 
relative velocities between particles are less than their average thermal velocities. More 
generally, some pairs of species require i's,t to be multiplied by a correct ing function of 



relative velocity, which reduces to unity in the linear limit we have taken (jSchunk and Nagy 



20091). 



We also work with the total momentum equation for the partially-ionized plasma as a 
whole: 

Dv ^ ^ 

p— = jxB-Vp + pG- (11) 



low-frequency Ampere's law: 



Faraday's law: 



J = -Vx5; (12) 
Po 



= -Vx^; (13) 



and generalized Ohm's law: 



E = -UixB + —j + —{jxB-Vpe], (14) 



- 15 - 



with 

<ye = — V (15) 

being the electron conductivity allowing for collisions with ions and neutrals (this is the 
same as the parallel conductivity in a single-fluid description of partially ionized plasma). 
Finally, the single fluid velocity is related to the neutral and ion velocities by 

PU = PnUn + PiUi, (16) 

where an electron contribution is neglected owing to small electron mass. 



3.4. One- Dimensional Linear Model 

For the present investigation, the model is reduced to one dimension. ID models are 
useful tools that provide insight by addressing the key physics of a problem in a simplifled 
way. Their results, and lessons from their development, can also pave the way for 3D 
modeling by establishing the most fundamental behavior of a system, proving techniques, 
and identifying potential problems before large resources are committed. If properly 
constructed, their results are often robust (to at least leading order) when compared with 
higher dimensional models. 

Our governing equations are simplifled by assuming horizontal invariance, linearizing 
the governing equations and neglecting the Hall term in Ohm's law (a valid approximation 
provided length scales remain longer than the ion skin depth, 6i = \/mi/ {poneC^), which 
is consistent with our results). Under these conditions, horizontal velocity perturbations 
decouple from vertical motions and do not compress the plasma or magnetic fleld, therefore 
taking the form of shear Alfven waves. One-dimensional modeling also requires that we 
take the equilibrium magnetic fleld direction as vertical, so that there is a single preferred 
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direction, and we neglect the effect of electron-neutral collisions in the neutral momentum 
equation in favor of ion-neutral collisions, owing to the much smaller momentum carried by 
electrons than ions. 

The linear assumption limits our model to early times during the response of the 
solar atmosphere to incident coronal waves, which is the focus of this study, and to wave 
amplitudes with 6B ^ Bq. Relaxing this assumption is an important target for future work 
motivated by this study and is discussed in Section [71 Other goals for future work include 
introducing second and third dimensions and more sophisticated modeling of the interaction 
betweens neutral particles and the plasma. These are also discussed in Section [71 

Using these simplifications, the equations governing shear Alfven waves become 

du Bo. 

at p 

= J^n,i [Ui - Un) , (18) 

dj fBo , 1 



dt dz^ \Ho fioae J 
There are computational advantages to using the three variables u, Un and j and these 
are described in Section [H The final step is to eliminate Uj, which is done with the aid of 
equation f[TB]l . Doing so, and making use of equation f lTU]) . the governing equations can be 
expressed as 



du. 



n 



{l^n,i + ^in){u-Un), (21) 



dt 

dj d"^ f Bo f On , A, 1 



dt dz^ \Po \ Pi J Poo-e J 

For the chromosphere, equations (l^U]) to (12^ are solved as written; for the corona, 
equation f[2T]) is dropped and the limits for fully ionized plasma at large magnetic Reynolds 
number (p„ — )■ and l/cTg — )■ 0) are applied to equations f l20|) and f l22|) before solving. 
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If the displacement current is retained in Ampere's law for the corona, then the 
right hand side of the ideal version of equation fl22]) is divided by a factor of 1 + c\/c^ 
(not shown above). Consequently, Alfven waves travel at the relativistic Alfven speed 
CA,rp\ = caI -\/1 + c \/c^ and therefore cannot exceed the speed of light (e.g. 



2008 



Russell 



Lysak and Song 



20101). In more extreme cases than are considered here, a relativistic treatment 
of the wave propagation may be required, however, for the semiempirical models studied 
in this paper, all Alfven speeds are nonrelativistic and the displacement current can be 
safely neglected (for reference, if ca = 0.1c, which is larger than the values we consider, the 
relativistic correction changes the wave speed by less than 0.5%). 

The goal of this study is to estimate transmission of wave energy to the chromosphere 
and quantify chromospheric damping, hence, the crucial aspects that must be treated 
realistically are wave speeds, electron resistivity and damping due to ion-neutral collisions. 
This leads to a conflict, because wave speeds and the effectiveness of ion-neutral collisions 
in wave damping both depend on the strength of the equilibrium magnetic field, which must 
be allowed to vary, even though we have neglected horizontal components when setting the 
equilibrium magnetic field direction to obtain equations (120|) to f l22|) . To get speeds and 
damping right, we therefore allow Bq to vary as a function of height when solving equations 
(1201) to (1221) . In consequence, the ID model is not self-consistent in terms of equilibrium 
field strength and direction, however we have been careful to evaluate the impact of this 
on our model, checking for artifacts, and we believe it does not unduly affects our results. 
For example, we have verified that: V • 5 = for perturbations; energy and momentum 
are both conserved and transported by appropriate vectors; and comparison to simulations 
using uniform field strength shows only the expected changes due to wave speed and 
damping rates. 
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3.5. Collision Frequencies 



Chromospheric ions are modeled as a mixture of protons, magnesium, silicon and 
iron, with the number density of heavy ions d etermined by quasi-neutrality and relative 



abundances taken from 



Kaye and Labyl (Il995l ). Heavy ions are assumed to be once-ionized 



where their numbers are significant (around the TMR). At chromospheric temperatures, the 
dominant interaction between neutral hydrogen and protons is resonant charge exchange, 
which transfers momentum with an effective collision frequency 



2.65 X IQ-^^T (1 - 0.083 log^g T)' Up. 



(23) 



Neutral hydrogen and heavy ions interact through non-resonant. Maxwell molecule 
collisions, giving 



2.11 X 10"^^ (1 + I/Am) 



-1/2 



(24) 



where Am is the atomic mass number of the ion. Combining these, the total neutral-ion 
collision frequency is 



(25) 



Note that the form of equation means i'H,i+ is insensiti ve to changes in th e comp osition 



of the heavy-ion fluid, 
with close agreement to 



he above frequencies a re given by 



Geiss and Burgil mm (GB). 



Schunk and Naevl mm (SN) 



A number of formulas for z/g- appear in the literature. Here, SN quote 



= 4.5 X IQ-^^T^/^ (1 - 1.35 X 10-^T) uh, 



(26) 



(11982 



(P)) quotes i'e-,H = 1-95 x 



GB quote u^-^h = 6.97 x IQ-^^T^-^uh and iPriestl 
]^g-i62ii/2^^ All three formulas give iye-,H > ^e-,i+ around the temperature minimum, 
meaning that electron-neutral collisions dominate electron resistivity at these heights, 
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where as electron-ion c ollisio ns dominate in the upper chromosphere and solar interior 
In the 



Fontenla et al. 



(120091 ) plage model, at heights where electron-neutral collisions are 
significant, the value for Ve-,H from GB is greater than the value from SN by a factor 
slightly less than 2, while the value obtained from P is less than the value from SN by a 
factor that can be as large as ten. We chose to use the SN formula because it produces 
the intermediate value where electron-neutral collisions are significant and is in reasonable 
agreement with the GB value. The standing uncertainty in the electron resistivity of the 
chromosphere's most neutral layers is noted, but fortunately has only a small impact on the 
conclusions of this work. 



4. Simulation Methods 

Equations fl20l) to fl22|) were solved by finite difference on a non-staggered grid, with 
the selection of primary variables (u, m„, j) avoiding grid decoupling. Chromospheric length 
scales are much shorter than coronal length scales (due to different Alfven speeds) so two 
different grid spacings were employed. Within each domain, d'^/dz^ was evaluated as a 
three-point centered finite difference. At the interface point, a five-point finite difference 
was used to maintain second order accuracy. 

The ion-neutral coupling time scale, l/(z/„ j + Vi^n)-, is on the order of milliseconds for 
the chromosphere and we will investigate waves with periods of order seconds, which makes 
equation (12T]) stiff. This is a significant consideration for selecting a numerical scheme and we 



have been successful using the odelSs MATLAB solver algorithm (iShampine and Reichelt 



19971 ). Accuracy is controlled by automatic adjustment of the time step to satisfy a relative 
tolerance of 10~^, an absolute toleranc e on velocities of 10~^° ms~^ and absolute tolerance 
on current density of 10~^° Am~^ (e.g. IPress et al.l 120071 ). 
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For this combination of finite differences and numerical scheme, numerical dissipation 
and dispersion become noticeable when about 30 points or less resolve one wavelength, the 
primary signature being an overshoot behind a wavepacket. Simulations presented in this 
paper have been checked for these artifacts, and comparison against simulations using fewer 
grid points shows them to be well converged. The two-grid approach was validated using 
simulations of Alfven waves in a fully ionized plasma, first in a uniform medium and then 
with a discontinuity in Alfven speed at the grid interface. In these tests, the change in grid 
spacing did not produce numerical reflections, provided the wave was spatially resolved in 
both domains. 

Simulations are driven by launching a downgoing wave in the corona, using a 
time-dependent boundary condition for j that allows upgoing waves to escape without 
reflection. Using the ideal limit of equation (122 p in the corona (p„ — t- and l/de — > 0) and 
splitting it into the sum of a contribution from the upgoing wave and specified downgoing 
wave: 

dj dji Bo d f duf 



dt dt ^ i^odz\dz ) ' ^^'^^ 

For ideal Alfven waves in a fully ionized corona, propagating in one direction, magnetic and 
velocity perturbations are related by b/Bo = ±u/va, where the plus sign is for downgoing 
waves and the minus sign is for upgoing waves. Thus, can be written in terms of total 
fields u and b, giving 



dj ^ dji 
dt dt 



1 






2 







dj_ 

dz 



(28) 



The spatial derivatives in equation ( l28l) apply to the upgoing wave and can therefore be 
evaluated as one-sided finite differences using points inside the domain, without causing 
numerical instability. We have used second order expressions and find that these produce a 
very clean outgoing boundary condition. 
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The incident wave is specified through dj^/dt in equation ( 128|) . For the simulations 
presented in this paper, is a single wavepacket that is both smooth and continuous: 
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(29) 



v 



where r is the driver duration and jo is the current-density amplitude of the downgoing 
wave. The resulting j{t) is initially zero, grows to j = jo at t = r/4, declines to change sign 
at t = r/2, reaches j = —jo at t = 3r/4 and returns to zero at t = r. 

Simulations were performed on four cores of an 15-750 processor (clock speed of 2.66 
GHz) with 4 GB of RAM. Run times ranged from less than one hour up to three days, the 
longer running simulations being those requiring the more grid points to resolve smaller 
length scales. 

5. Case Study: Plage Atmosphere, 1 Second Driver 

A case study is now presented, using a plage model and a 1 second duration driver. 
This allows for detailed discussion of the evolution, the features of which are common 
to all drivers and atmospheric models. Wave transmission to the chromosphere and 
damping within it are also quantified. Having established key behaviors in this section, a 
parameter study will be presented in Section [6] which examines the effects of wave period 
and atmospheric model. 

We begin by describing the how the wave evolves in the case study. First, the driver 
generates a downgoing Alfven wave in the corona that partially reflects from the transition 
region. The majority of wave energy (78.0%) is reflected to form an upgoing Alfven wave 
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in the corona, which escapes from the simulation, with the remaining 22.0% transmitted to 
the chromosphere. 

Once in the chromosphere, the transmitted wave evolves as shown in Figure 121 which 
plots wave energy density (the sum of magnetic and kinetic energy densities) as a function of 
height at 2.5 second intervals. The wave propagates downwards through the chromosphere 
with very little additional reflection because the width of the pulse (approximately 100 km 
at its widest; t = 2.5 s) is less than the length scale of Alfven speed gradients. Here, 

— 1/2 

ca ~ Bufj changes primarily through hh, which has a scale height of about 230 km in 
the upper chromosphere (see equation ([2]) and Figured]), thus the Alfven speed length 
scale would be approximately 460 km if magnetic field were ignored. In fact, magnetic 
field strength increases with decreasing altitude, partly compensating for the decrease in ca 
caused by increasing nn, so the true length scale of Alfven speed is longer than 460 km. It 
follows that the spatial scales of the wave perturbations and equilibrium are well-separated 
in this case, hence, the chromospheric wave evolves as in a WKB description and is not 
significantly reflected by Alfven speed gradients. 

Decreasing Alfven speed does affect the Alfven wave by slowing the pulse and 
altering its shape. The trailing edge travels faster than the leading edge, thus wavelength 
decreases with height, which in turn increases maximum wave energy density (although 
height-integrated total energy density is conserved). 

Of the wave fields, current density is most affected by the shortening wavelength. 
Magnetic-field perturbations increase as the wavelength shortens (to maintain height- 
integrated magnetic energy if damping rates are not large), and current density is 
proportional to magnetic-field perturbation divided by wavelength. Thus, current density 
would be inversely proportional to the square of wavelength in the absence of damping. In 
this simulation, which includes physical damping, maximum current density at t = 50 s is 
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28.9 times greater than maximum current density at t = 2.5 s. 

Figure [3] plots energy densities integrated in height over the chromospheric domain 
and expressed as a fraction of the energy injected to the simulation by the incident coronal 
Alfven wave. One second into the simulation (equal to the travel time from the simulation's 
upper boundary to the transition region) the energy in the chromosphere sharply increases 
as the wave enters the chromosphere. In this case, 22.0% of incident wave energy is 
transmitted to the chromosphere (evaluated using energy at the end of the simulation). 
This is high enough to be of considerable interest. 

As expected from analysis of Figure [21 Figure [3] shows that very little energy is returned 
to the corona by subsequent reflections inside the chromosphere, the total energy in the 
chromosphere remaining very nearly constant after the initial rise. 

From the first moment when the Alfven wave is in the chromosphere, it loses energy. 
The decrease in wave energy is gentle when the wave is in the upper chromosphere, then 
rapid as it passes through the TMR, and more again gentle below this. Energy lost from 
the Alfven wave becomes heat, and 54.6% of (non-equilibrium) energy in the chromosphere 
at the end of the simulation is heat that has been produced by wave damping. 

To identify the source of this damping, we note that heating occurs at a rate 



where the first term on the right-hand-side corresponds to ion-neutral friction and the 
second term is Ohmic heating due to electron resistivity. Figure |3] plots both heat 
contributions. When the distinction is made, it is clear that wave damping is primarily due 
to ion-neutral friction, although electron resistivity plays a non-negligible role towards the 
end of the simulation. 

Figure H] examines the locations of heating, with Figure H] (top) showing the thermal 




(30) 
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energy produced by ion- neutral friction and electron resistivity as a function of height. It 
confirms that ion-neutral friction is the dominant heating mechanism throughout most 
of the chromosphere, with electron resistivity taking over below about 0.4 Mm. Very 
noticeably, the thermal energy deposited by wave damping peaks dramatically between 
about 0.4 to 0.9 Mm, which corresponds to the TMR. Further investigation has shown that 
ion-neutral friction peaks in this layer because it is where the ion-neutral coupling time 
is longest (i.e. i'H,i+ is smallest). This means that the relative velocity between ions and 
neutrals is greatest here, and this increases the heating rate in equation fl30l) . (The increase 
in relative velocity enters equation flSU]) as a square and therefore more than compensates 
for the reduction in VH,i+-) There is also a plateau of heating in the upper chromosphere, 
again due to ion-neutral friction. Meanwhile, below the TMR, electron resistivity becomes 
significant, partly due to the many-fold increase in j that we have described. 

Figure H] (bottom) shows the equivalent temperature increase that corresponds to the 
heat deposited by this wave (energy deposited per particle). We do not attempt to address 
the thermodynamics and the loss terms that would ultimately govern the temperature 
profile, however it is instructive to make a simple direct estimate which serves as a crude 
indication of the likely locations and extent of temperature changes. Although considerably 
more thermal energy is deposited at the TMR than in the upper chromosphere, heating the 
upper chromosphere produces a larger change in temperature because of the much lower 
particle density. Indeed, for a modest amplitude coronal wave with a magnetic perturbation 
of 5% of the background field {5B/Bq = 0.05), a single 1 second pulse heats the upper 
chromosphere by a peak amount equivalent to a temperature increase of 137 K. A secondary 
peak is also present close to the temperature minimum, where the heating is equivalent to 
11.4 K in spite of the relatively high densities. Values for other wave amplitudes can be 
obtained by scaling energy density and temperature as {5B/ BqY since we are using a linear 
model. These temperature increases are large enough to be interesting for chromospheric 
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emission, discussed further in Section [71 



6. Parameter Study 

Having established that transmission of wave energy from the corona to the 
chromosphere can be significant and that chromospheric waves can be substantially damped 
(heating the chromosphere) we now investigate sensitivity to different wave periods and 
atmospheric models. 



Chromospheric models were based on the semiempirical models of iFontenla et al. 



(120091 ). with magnetic field strength models provided by equation (jS]). The six models are 
therefore determined by the parameters shown in Table [TJ Simulations were performed 
for each of these models, with wave periods between 0.5 and 30 seconds. In all cases, 
simulations were qualitatively similar to the case study but differed in the fraction of energy 
transmitted to the chromosphere and the effectiveness of wave damping. 

Table 1. Chromospheric models used for parameter study. 



Description 


Fontenla Model ID 


Bphot (T) 


Bcor (T) 


Umbra 


1006 


0.2 


0.05 


Weaker Umbra 


1006 


0.2 


0.025 


Penumbra 


1007 


0.1 


0.02 


Plage 


1004 


0.05 


0.01 


Facula 


1005 


0.05 


0.01 


Bright Facula 


1008 


0.05 


0.01 
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Height (Mm) 

Fig. 2. — Snapshots showing evolution of energy density as a function of height above the 
photosphere for case study of a 1 second period wave in the plage model. See the electronic 
edition of the Journal for a color version of this figure. 




Fig. 3. — Energetics of the chromospheric domain for case study of a 1 second period wave in 
the plage modeL The wave is launched in the corona at t = s and enters the chromosphere 
at if: = 1 s. See the electronic edition of the Journal for a color version of this figure. 
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0.09 



■ Ion-Neutral Friction 
Electron Resistivity 




1 1.2 
Height (Mm) 



1.8 



Fig. 4. — (top) Thermal energy deposited by a 1 second period wave in the plage model, 
through ion-neutral friction and electron resistivity, (bottom) Equivalent temperature in- 
crease produced by damping of the wave packet in the chromosphere (thermal energy per 
particle). In (top) and (bottom), values correspond to a wave with coronal amplitude 
5B/Bq = 0.05 and scale with (SB/Bq)"^. See the electronic edition of the Journal for a 
color version of this figure. 
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6.1. Energy Transmission Prom Corona to Chromosphere 

Figure [5] shows energy transmission to the chromosphere for each simulation, evaluated 
as the energy that is located in the chromosphere at late times divided by the energy 
injected to the simulation by the incident coronal Alfven wave. 

The most striking feature is a strong dependence on wave period, transmission being 
greater for shorter period (higher frequency) waves. Looking at the plage model, energy 
transmission is 24.1% for a 0.5 s period wave, 22.0% for a 1 s period, 13.9% for a 10 s 
period and 7.6% for a 30 s period. This variation occurs because wave reflection is not 
determined solely by transmission across the transition region, but also by wave properties 
in the chromosphere. 

The transmission analysis of Section |2] is a good approximation when wavelengths 
are much shorter than coronal and chromospheric speed gradient length scales. In this 
short-period limit, waves propagate without reflection once in the chromosphere, as was 
seen in the 1 s period case study. This places an upper limit on the energy transmission 
which, using equation ([7]), is 23.2% for the plage model (1 MK corona and 6700 K 
neutral-dominated chromosphere). In fact, the shortest period waves can beat this limit, 
which we attribute to the presence of a thin layer of fully-ionized plasma at the base of 
the transition region, equation ([7]) having shown that transmission to an ionized plasma is 
higher than transmission to a neutrally-dominated plasma. 

The assumption of uniform media breaks down when waves in the upper chromosphere 
have a wavelength comparable to or larger than the speed gradient length scale, in which 
case waves become subject to additional partial reflections that occur continuously. In the 
plage model, 1 second period waves have a wavelength that is indeed much smaller than 
chromospheric scale heights, so transmission approaches the limit described by equation 
(I7j). Longer period waves, however, are reflected in the upper chromosphere, and the longest 
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period waves are subject to the greatest additional reflection, as evident in Figure |5l 

Since wave speed decreases with height (as a result of rapidly increasing density), 
wavelengths decrease significantly as waves propagate through the chromosphere, as seen in 
Section [51 Scale heights do not change as significantly, e.g. staying in the range 120 km 
(lower chromosphere) to 230 km (upper chromosphere) in the plage model. Therefore, waves 
cease to reflect as they reach lower altitudes, hence reflection within the chromosphere only 
happens in its upper reaches. 

Upgoing waves produced by continuous reflection in the upper chromosphere partially 
reflect at the transition region to form new downgoing waves. Therefore, wave energy is 
not easily lost from the chromosphere and multiple reflections can occur. These partial 
reflections act to extend the downgoing wave by superposition behind the primary pulse. 

Transmission also varies for different solar features. The umbral models are notably 
different because they have a cooler chromosphere (minimum temperature of 3590 K) and 
are therefore denser. They also have a lower transition region than other models, at about 
1 Mm altitude, and do not have the upper temperature plateau found in other models (e.g. 
see Figured] (middle)). This reduces transmission in two ways. 

First, temperature contrast between corona and chromosphere - and therefore density 
contrast - is greater for umbral models than for models of other solar features. The 
basic analysis of Section |2] suggests that energy transmission across the transition region 
scales as the square root of the temperature just below the transition region. Given that 
temperatures below the transition region are about 4400 K in the umbral model and 
6700 K in the plage model, we therefore expect energy transmission for umbral models to 
be approximately 0.8 times that for the plage case (for short period waves and at common 
coronal temperatures). This is in line with the reduction seen in simulations using a 1 s 
period driver, from 22.0% transmission for plage to 15.4% transmission for umbra and 
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18.0% transmission for the weaker umbra model. 

The second consideration is that since the hydrostatic scale height is proportional to 
temperature, the cool umbral chromosphere has a much shorter density scale height (about 
88 km) than do other solar features (e.g. 230 km in the upper plage chromosphere). Waves 
are therefore subject to greater reflection within the umbral chromosphere than within any 
other feature. They are affected for longer durations, and reflection affects shorter period 
waves, because a shorter wavelength is needed to turn off reflection. It follows that the 
greatest differences in transmission to the chromosphere, when comparing different solar 
features, occur at longer periods (e.g. for a 30 s period wave, transmission is 7.6% in the 
plage model and a meager 1.7% in the umbra model). 

The differences between the umbra and weaker umbra models are due to the magnetic 
field model. Both have a field strength of 0.2 T at Mm, but 0.025 T is applied at the 
top of the weaker umbra chromosphere, compared to 0.05 T at the top of the umbra 
chromosphere. Thus, magnetic field strength changes more over the same distance and 
therefore has a shorter length scale in the weaker field model (the factor is 1.5). Looking 
at equation ([2]), increases in field strength with depth counteract decreases in Alfven speed 
caused by increasing density, so Alfven speed decreases more gradually in the weaker umbra 
model than in the umbra model. Thus, in the weaker field model, waves experience less 
reflection within the chromosphere than for the stronger field model, increasing overall 
transmission. 



6.2. Chromospheric Damping 

Chromospheric damping for different drivers and atmospheric models is shown in 
Figure O which plots the percentage of energy in the chromosphere that is heat, at the end 
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of each simulation. Contributions from ion-neutral friction and electron resistivity are both 
identified. 

For 1 s period waves, between 36.6% and 79.2% of wave energy injected into the 
chromosphere has damped by the end of each simulation. This damping is primarily due to 
ion-neutral friction, with the greatest losses occurring around the TMR, although electron 
resistivity is responsible for up to a quarter of the total damping of 1 s waves in an umbra. 

Damping is very sensitive to wave period. In the plage model, 76.8% of wave energy 
entering the plage chromosphere is damped for 0.5 s periods; at 1 s periods damping falls 
to 54.6%; and at 10 s periods, only 1.8% is damped. Thus, waves with periods of about 
10 seconds or greater pass through the chromosphere to the interior relatively undamped, 
while waves with periods of about 1 second or shorter are strongly damped and heat the 
chromosphere effectively. 

A simple estimate of the damping time due to ion-neutral collisions illuminates the 
sensitivity of damping to wave period. Using equations ( TT8l) and ( l30l) . the wave loses energy 
through ion-neutral friction at a rate 

dE^nave. , \2 /I dUn\ /oi\ 



at \Vni at 



The derivative dun/dt can be approximated as UnlT^ave where r^ave is the wave period. Also 
note that kinetic energy of neutrals makes up almost all kinetic energy in the chromosphere 
and that kinetic and magnetic energies are approximately equal once waves are free to 
propagate downwards (confirmed against our simulations), thus, pnU^ ~ E^ave^ and hence, 

~ (32) 

making the energy damping time 

'^damping ~ ^n,iTwave- (33) 
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Although this estimate is missing a coefficient of order unity (which will depend on the 
shape of the wave), we see that ion-neutral friction is most effective for short period waves 
and small values of z/„ j (longer coupling times). Of these two parameters, damping is more 
sensitive to wave period, which enters as a square. 

Comparing different solar features across Figure El damping is strongest in the umbral 
models, where almost 80% of wave energy entering the chromosphere is converted to 
heat. Ion- neutral damping is more effective in umbras than for other features because 
the neutral-ion coupling time is longer here - typically one order of magnitude longer 
when compared with a plage model. This allows larger relative velocities between ions 
and neutrals, which results in greater heating. Damping due to electron resistivity is also 
substantially stronger in umbras than other solar features. Comparing the umbra and plage 
models, electron resistivity (l/cTg) at 0.75 Mm altitude is similar in both models, however 
it sharply diverges below this height, becoming as much as 250 times greater in the umbra 
model and this causes greater damping of chromospheric waves (see equation (|30|) ). 



7. Conclusions 

Our studies have shown: 

1. Coronal Alfvenic waves can transmit a significant fraction of their energy to the 
chromosphere. Energy transmission is greatest for short period waves and can 
exceed 20% per incidence for waves with second or subsecond periods. Cumulative 
transmission will be higher if coronal waves are trapped in closed coronal structures 
(e.g. closed magnetic field), leading to multiple incidences. 

2. Waves in the chromosphere are subject to damping by ion-neutral friction and 
enhanced electron resistivity, with ion-neutral friction dominating in most cases. 
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30% 




Driver Period 



Fig. 5. — Energy transmission coefficients for waves incident from the corona. For each 
atmospheric model, transmission is shown for driver periods spanning 0.5 to 30 seconds. See 
the electronic edition of the Journal for a color version of this figure. 
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Fig. 6. — Percentage of chromospheric wave energy converted to heat by ion-neutral friction 
and (neutral enhanced) electron resistivity. For each atmospheric model, damping is shown 
for driver periods spanning 0.5 to 15 seconds. See the electronic edition of the Journal for a 
color version of this figure. 
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Damping is highly sensitive to period: waves with subsecond periods are damped 
effectively (for 0.5 s periods up to 91% of wave energy entering the chromosphere is 
converted to heat), whereas waves with periods longer than 10 s pass through the 
chromosphere with relatively httle damping (less than 7.2%). 

3. Chromospheric damping of downgoing waves delivers most energy to the temperature 
minimum region (TMR), heating it directly. Energy is also deposited throughout the 
upper chromosphere, where temperature increases are expected to be larger than at 
the TMR because number densities are much lower. 

These findings suggest that flare-generated MHD waves may be a viable means 
of transporting flare energy from the corona to the chromosphere and heating the 
chromosphere. They would have greatest effect if a substantial part of wave power were at 
second or subsecond periods. Periods typically associated with solar flares, that one might 
expect to flnd in flare-generated waves, range from subpulses in X-ray, radio and visible 
light observations (order of 100 ms), through elementary flare bursts (tens of seconds), to 
the duration of the rise phase (several minutes). There is, we conclude, reasonable overlap 
between observed flare time scales and the wave periods that have potential to heat the 
chromosphere. 

For a case study of a single 1 second wave pulse in a plage model, a modest magnetic 
fleld perturbation of 5 G on a 100 G coronal background fleld produced temperature 
increases of up to 137 K in the upper chromosphere and 11.4 K at the TMR. Heating 
could be considerably amplified by increases in perturbation amplitude or background field 
strength. If a train of such pulses were sustained for a minute, the resulting temperature 
changes would certainly be significant. 

It is interesting that the bulk of thermal energy is delivered to the TMR (although 
temperature rises there are less than those in the less dense upper chromosphere). This 
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profile of heating is very different to that for electron beams, which lose their energy higher 
in the chromosphere. Waves therefore offers a direct mechanism for heating the TMR, 
where the bulk of white light flare emission is formed, that should be considered as a 
complement to backwarming. 

Wave energy in the chromosphere that is not lost to heat is available to particle 
acceleration, since Alfven waves with fields that vary perpendicular to the background 
magnet ic field can carry magn etic-field- aligned electric fields created by inertial or kinetic 



effects ( jStasiewicz et al. 



2OOOI ) or electron resistivity. They could, therefore, potentially 
contribute to producing fluxes of energetic electrons in the tightly focused X-ray footpoints 
of a solar flare. This may help to relieve theoretical difficulties associated with electron- 
beam-only coUisional-thick-target models, such as apparent requirements for large particle 
numbers, strong anisqtropy and high beam densities in the corona that would be unstable 



( iFletcher and Hudson 



2008 



Brown et al 



2009 



Krucker et al 



20111 ) 



Electron acceleration can only take place in the chromosphere if there is a population 
of electrons for which energy losses to collisions can be overcome by energy gains between 
collisions. This requires some electrons to have speeds significantly greater than typical 
quiet chromosphere thermal speeds (for which collisional losses ar e severe). Two ways in 
which this can be achieved are by: (i) heating the chromosphere (iFletcher and Hudson 
20081 ). e.g. by damping of flare-generated MHD waves; or, (ii) by having electrons present 
that previously precipitated from the corona and ret ain sufficient veloci ty following each 



collision to be reaccelerated before the next collision f lBrown et al 



20091). 



In the case of directly accelerated chromospheric electrons, one can imagine a two-stage 
evolution, in which the upper chromosphere is initially heated by waves, then energy is 
increasingly given over to acceleration of chromospheric electrons as the high energy part of 
the electron distribution becomes fast enough for acceleration to overcome collisional losses. 
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A similar picture was presented in the flare burst model of iLiu and Fletcherl ( 120091 ) . This 
would be consistent with observations that some flares show increased thermal emission up 
to minutes before the rise phase of hard X-ray (HXR) emission. 

In such scenarios, heating and electron acceleration would draw energy from the same 
source (waves), giving rise to a relationship between heating rate (hence time derivative of 
thermal emissio n) and HXRs. Such a relationship may produce correlations similar to the 



Neupert effect (INeupert 



1968h. which i s obse rved in a significan t fract ion of flares (about 



half those studied by iMcTiernan et al. 



Veronig et al. 



( 120021 ) )■ There is, however. 



mm and 

an advantage to wave models, that if energy were delivered to the chromosphere (at least 
in part) by waves, then thermal and HXR emission would not need be as closely related 
in all cases as is required by electron-beam-only models; subsequently, flare models that 
include wave fluxes may be able to better explain non-Neupert flares. Unfortunately, the 
exact form of relationship that would arise from a wave model, and the circumstances under 
which it would hold are not known at this time. 

A further point to discuss is that not all wave energy is lost before downgoing waves 
pass into the solar interior. This provides a mechanism by which energy and momentum 
can be transferred from the corona to the i nterior during a solar flare , as required to exc ite 



1998 



Kosovichev 



sunqu akes (flare-generated seismic waves) (iKosovichev and Zharkova 
20061 ). Some aspects of su n quake excitation by direct wave coupling were previously 



presented by 



Hudson et al 



(I2OI2I ) and, although such a connection remains speculative, 
the present work shows waves with periods of several seconds have the best transmission 
characteristics to enter the chromosphere and pass undamped into the solar interior. 

This study has used linear models based on the preflare atmosphere to establish initial 
transmission and damping coefficients, subject to a variety of simplifying assumptions that 
limit the validity of these models. These have established that MHD waves are of interest to 
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the flare chromosphere and our results justify further investigation using more sophisticated 
models. Limitations of the present models and their possible consequences on the results 
are now discussed. 

The initial transmission and damping of waves with periods of 1 second or less is so 
efficient that it will quickly lead to expansion of the heated layers, changing densities, 
increased radiative emissions and changes to the ionization state; none of the above having 
been considered by the present investigation. Thus, non-linear studies are required to 
properly determine longer term evolution, even over the few minutes of the flare impulsive 
phase. Such studies must also give careful consideration to conduction from the heated 
upper atmosphere to the TMR, and optically thick radiative transfer and losses (which 
are sensitive to spectral lines), effects that were avoided in the present study by the linear 
assumption. 

During longer term evolution, we expect that chromospheric evaporation will decrease 
the density ratio between the corona and chromosphere, increasing transmission to the 
chromosphere beyond that found in this study. Thus, waves may become an increasingly 
effective mode of energy transport as the atmosphere evolves. At the same time, changes 
in density and temperature will alter collision frequencies. In the TMR, heating will 
lead to increased ionization, increasing z/„^j, and the proton contribution to Un^i will also 
increase through its temperature dependence (see equations fl23|) and ( 124|) ). This will 
act to decrease the effectiveness of wave damping by ion-neutral friction. On the other 
hand, density decreases due to expansion of the heated layers will act to decrease j 
and therefore increase the effectiveness of wave damping. Electron resistivity will also be 
affected: increased ionization will increase the ion contribution and decrease the neutral 
contribution, while direct temperature increases and density reductions decrease both. The 
full picture is therefore likely to be a series of positive and negative feedbacks that increase 
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wave transmission but ultimately quench wave heating. 

The present study is additionally limited by a ID assumption. It is desirable that 
this should be relaxed to consider incident waves with finite transverse wavelengths (which 
must be present given that the corona is highly structured) and magnetic field geometry 
consistent with height variation of the magnetic field strength (as discussed in Section 



13. 4[ the present model allows the equilibrium magnetic field strength to vary with height, 
although horizontal components of the equilibrium field are neglected; this does not appear 
to introduce undue artifacts, however, it means that the model is not fully self-consistent). 
Non-linear wave amplitudes should also be considered for fiares, leading to wave steepening, 
with increased dissipation, altered transmission properties, and generation of turbulence. 
The last area we identify for future improvement is modeling of the different particle 
species and their collisions. Ion-neutral collisions were assumed elastic in this work, which 
is reasonable when the relative ion-neutral fluid velocity is less than typical thermal speeds. 
It is possible, however, that studies of the longer term atmospheric response to wave 
energy will show this con dition bre aks down, in which case coUisional or Alfven ionization 



may play some role (e.g. 



Lai 



2001 



and references therein). Furthermore, in a rapidly 
evolving atmosphere, populations ought to be solved explicitly in time without equilibrium 
assumptions. 

In conclusion, this work has shown that flare-generated MHD waves can transport 
energy from the corona to the chromosphere and heat the chromosphere, with short period 
waves (a few seconds or less) being most effective. They are particularly interesting as 
a mechanism for heating the TMR, thereby contributing to white light emissions. This 
paves the way for more sophisticated modeling that may establish a viable counterpart to 
electron-beam flare models. 
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